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' Abstract. We describe some semiclassical spectral properties of Harper-like operators, i. e. 

of one-dimensional quantum Hamiltonians periodic in both momentum and position. The 
spectral region corresponding to the separatrices of the classical Hamiltonian is studied for 
the case of integer flux. We derive asymptotic formula for the dispersion relations, the width 
' of bands and gaps, and show how geometric characteristics and the absence of symmetries of 

the Hamiltonian influences the form of the energy bands. 

' 1. Introduction 

o\ . 

In the present work we are going to describe certain asymptotic spectral properties of the Harper-like operators. 
Such operators appear as follows. Let H(p,x) be a real-valued real-analytic one-dimensional classical Hamiltonian 
periodic in both momentum and position: 

H(p + 2n,x) = H(p : x + 2n) = H(p,x), p,x e R. 

O 

will be called the Harper-like operator associated with H. In one of the simplest cases, when H(p,x) = 2 cos p + 
2a cosx, one has ifi^f) (x) = f(x +h)+ f(x — h) + 2acosx f(x), i. e. Hi, is the Harper operator on the real line. 

Such operators appear in the study of the motion of electrons in a two-dimensional periodic potential subjected 
to a uniform perpendicular magnetic field [1,2], and the parameter h may be expressed through the parameters of 
the system differently in the cases of strong and weak magnetic fields [3]. The periodic magnetic systems show a 
rich spectral structure depending on h, which in general expresses a relationship between the magnetic and electric 
fields; for example, in the strong magnetic field one can put h = (Im/L) 2 , where Im is the magnetic length, and L is 
the characteristic size of the period lattice [4]. The Hamiltonian Hf, commutes with operators Tj, j — 1,2, defined 
by 

(7i/)(x) = r^/W, (T 2 f)(x) =f{x + 2n)- 
they obey the equality T\T2 = exp(47T 2 / / \i)TjT\ and commute iff r\ := 2%jh G Z; this number tj is usually referred 
to as the number of magnetic flux quanta through an elementary cell or simply flux. In the case of integer 77 
one can apply the usual Bloch theory and show that the spectrum consists of Tj bands [5], so that each band 
is the value set of the corresponding dispersion relation E(k,h), where k = (&i,&2) £ [—n,n) x [—7Z,n) is the 
(vector) quasimomentum, and for any k there exists a (generalized) eigenfunction ^(x, k, h) of H/„ Hf^V{x : k, h) = 
E(k,h) l ¥(x,k,h), satisfying the Bloch-periodicity conditions 

Tj¥(x,k,h)=e? k W(x,k,h), 7 = 1,2. (1) 

The situation with non-integer but rational 77 = N/M (N and M are mutually prime integers) can be reduced to the 
previous case by enlarging the unit cell of the period lattice (27r,27r) — > (2nM 1 2n). For irrational tj the spectrum 
is rather complicated and can include parts of Cantor structure as was predicted in [2, 6] and then analytically 
justified for the Harper operator in [7, 8]. 

We will study the asymptotics of the spectrum for ///, for large integer 77, which corresponds to the semiclassical 
limit (h — > 0). Spectral bands in different parts of the spectrum have then different asymptotic behavior [5, 1 1]. To 
illustrate this, we consider first the Harper operator. To be more definite, assume < a < 1, then, for any 8 > 0, 
the bands lying below (— 1 + a — 8) and over (1 — a + 8) correspond to finite classical motion and have the width 
o(h°°), while the bands inside the segment [—1 + a + 8,1 — a — 8] correspond to open classical trajectories and 
have the width ~ h with o(/i°°)-gaps between them [9]. These estimates are non-uniform with respect to 8, and they 
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Figure 1 . Space of trajectories as the Reeb graph, (a) Level curves of the Hamiltonian H in 
the unit cell, (b) The corresponding Reeb graph. Separatrices are bold and correspond to the 
branching points of the Reeb graph. 

do not describe the bandwidth asymptotics near the critical points ±(1 — a), where the eigenfunctions undergo a 
transition from localized to extended behavior [10]. For Hamiltonians of a more general form the spectral structure 
is suitably described by the Bohr-Sommerfeld quantization rule and can be illustrated with the help of the so-called 
Reeb graph technique [5, 1 1], so that critical values of H divide the spectrum into parts with different asymptotic 
behavior (we illustrate this in section|2j, and the transitions between these parts semiclassically correspond to the 
separatrices of the classical Hamiltonian. 

This transient region was earlier studied from different point of view [7, 8, 12] for the Harper operator and some 
closed to it Hamiltonians with high-order symmetry. From the other side, the absence of symmetries appears to 
be the generic situation, which influences the structure of spectral gaps and dispersion relations, which are very 
different from the ideal model of the Harper operator; this leads to new phenomena in the transport properties of 
the periodic magnetic systems, see the review [15]. (This covers not only the Harper-like operators, but also the 
periodic Schrodinger operators, where symmetry properties are in connection with the absolute continuity of the 
spectrum [16].) Our aim here is to provide a uniform approach to studying the band- and gapwidth transition and 
construction of the dispersion relations. In particular, we are interested in the way how the shape of the classical 
trajectories influences the shape of the energy bands. We study in details only some basic cases (sections|4]and[5}, 
but even these simple examples show a rich structure of the dispersion relations; in particular, it turns out that the 
shape of the energy bands as well as the ratio bandwidth/gapwidth in non-symmetric cases is very sensitive to the 
parameters of the system; this is discussed in section[6] 

It is worthwhile to emphasize that even the spectral problem for a single two-dimensional periodic system can 
lead to a number of Harper-like operators with essentially different spectral properties. We illustrate this with the 
example of the two-dimensional Landau operator with a periodic electric potential v; in the case of strong magnetic 
field the Schrodinger operator is 

1 / d \2 h 2 d 2 

L '= o\ ~ ih ~r~ +X2 ) +ev(x 1 ,x 2 ), (2) 

2 V dx\ ) 2 dx 2 

with small h and £ [4]. The spectral problem for L is reduced to a series of spectral problems for the Harper-like 
operators associated with the classical Hamiltonians 

L n {Y u Y 2 ) = {n + ^)h + eJ (^~(2n + l)hA Y )v(Y u Y 2 ) + O{e 2 ), neZ+, (3) 

where Jq is the Bessel function of order zero, = d 2 jdY 2 + d 2 jdY 2 (we explain this reduction in Appendix), 
and (Y\ , Y 2 ) are canonically conjugate. The Hamiltonians L„ have, generally speaking, trajectories of different kind 
for different n, which results in a difference between their spectra. 

2. Regular Bohr-Sommerfeld rules and the asymptotics of the dispersion relations 

In this section, we recall some simple constructions which are useful for estimating the dispersion relations [5, 
11]. 

2.1. Reeb graph. Throughout the paper we assume that all the critical points of H are non-degenerate. Define 
on the torus T 2 X := M 2 x /(2n,2n), which will be called the reduced phase space, an equivalence relation ~ by 
x\ ~ x 2 {x\ and x 2 lie in a connected component of a level set of H}, then the set G :— T 2 }x / ~ is a certain 
finite graph called the Reeb graph of H. The end points of the graph correspond to extremum points of the 



3 



E(k.,k,,h) 




(a) 



(b) 



Figure 2. Spectrum obtained with the help of the regular Bohr-Sommerfeld rules, (a) Point of 
the Reeb graph, which correspond to the trajectories satisfying the quantization condition; (b) 
Semiclassical dispersion relations outside the separatrix region 



Hamiltonian while the branching points correspond to saddle points and separatrices (see illustration in figureQ. 
It is natural to distinguish between edges corresponding to open trajectories on R 2 X and, respectively, to non- 
contractible trajectories on T^. (these edges of the Reeb graph will be referred to as edges of infinite motion) and 
to closed ones on the plane and contractible ones on the torus (edges of finite motion). This graph provides a good 
illustration for the structure of the trajectory space of H as well as for the semiclassical spectral asymptotics. 

2.2. Finite motion. Consider a spectral interval A = [E 1 ,E"] such that for each EGA the level set H = E on the 
reduced phase space T 2 X consists of a single closed trajectory. In other word, the corresponding part of the Reeb 
graph must be an edge of finite motion; in the example shown in figure ^ the interval [H + + 8, maxH], 8 > 0, 
satisfies this condition; the corresponding edge is e\ . The regular Bohr-Sommerfeld rules select from the family of 
all these closed trajectories a discrete family of trajectories satisfying the quantization condition 

1 / 1 

<p pdx € Z. 

2nh J 1 2 

(The construction of these trajectories can be interpreted as a selection of a certain discrete subset on all the edges 
of finite motion of the Reeb graph, see figure^). Each of these trajectories, A, implies a quasimode (E, E € K, 
Y(x,h) GL^(K), with \jf microlocally supported by A, E — H\ A + 0(h 2 ), and \\(H h — £)i^||/|| vl — o(h°°), which 
means that dist(spec Hi u E) = o{h°°). Clearly, to each of such E there corresponds a whole family of quasimodes. 
Namely, if A is a closed trajectory satisfying the quantization condition and (E, is the corresponding quasimode, 
then the trajectory Aj : = A+ (2k j\,2nj2), j = (71,72) G also satisfies the quantization condition and produces 
the quasimode (E,\j/j) with the same value E and \j/j(x,h) = exp(2nij ix / h)y/(x ~ 2nj2, h) = T x n T 2 12 y(x,h). Let 
us try to satisfy the Bloch conditions by a function of the form 

¥(x,k,h)= £ Cj(M)v/jM)= E Cj(k,h)T- h T- j2 V (x,h). (4) 

}=UiJ2)ez 2 j=(;i,; 2 )GZ 2 

Clearly, the coefficients Cj must solve a linear system: C ;i +ij 2 = Cj lt j 2 e l , Cj u j 1+ \ ~ Cj l j 2 e' k2 , 71,72 € Z, there- 
fore, Cj = ce'^ for some constant c, and the corresponding Bloch quasimodes have the form 

¥(x,k,h) = c £ e i ^e 2Kij ^ ll \i/(x~2Kj 2 ,h). (5) 

j=(ilJ2)GZ 2 

Therefore, to construct a Bloch quasimode in the form Q we do not need to satisfy any additional conditions about 
relationship between E and k, which results in semiclassically constant dispersion relations. 

If for each £eA the level set H contains several closed trajectories (and intersect several edges of finite motion 
(for example, it is the interval [min//,//_ — 8] in the example of figure^, the procedure described above is still 
applicable and gives the asymptotics of the spectrum up to o(h°°), but, in the case of some symmetries between 
families of trajectories degeneracies of the eigenvalues may occur; computation of their splitting is much more 
delicate [5, 17]. 

2.3. Infinite motion. Now consider the spectrum in the interval A = [E',E"} assuming that the correspond- 
ing region of the Reeb graph consists of two edges of infinite motion (the interval [H- + 8,H + — 8] and the 
edges <?2 and <?3 in the example of figure [0, i.e. for any E 6 A the level set H = E on T 2 consists of two 
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non-contractible trajectories, and in R 2 there are two families of open periodic trajectories. Clearly, there ex- 
ists a vector d = [di^do) G Z 2 , non-divisible by any other vector with integer- values components, such that 
for any of these trajectories, A = yp = P(t), x — X(f)), there exists a nonzero number T = T(A) satisfying 
(P,X)(t + T) = (P,X)(t) + 2nd for all t G R. To simplify the calculation we assume that d = (1,0) (otherwise 
one can choose new canonical coordinates on the phase plane, such that d = (1,0) in these new coordinates). 
Obviously, there are two families of trajectories corresponding to different edges of the Reeb graph: for the first 
family, the number T can be chosen positive, while for the another one it must be negative; we denote these families 
by e + and e~ respectively. 

Each trajectory A* = (p =P(t),x = X(t)) G e ± implies a quasimode (E ± ,\l/ ± ), E ± = H\ A ± +0(h 2 ), satisfying 
T 2 y/ ± (x,h) = \jf ± (x + 2n,h) =exp(/5 ± (A ± )//i)v/ ± (x,/i), where 

S ± (A ± )= f T p ± (t)dX ± (t). (6) 
Jo 

Clearly, the correspondence = H\ A ± «-* A* is one-to-one. Therefore, 5 ± can be considered as a function of 
E , S = S ± (E ± ). As this dependence is continuous and monotonic, it can be inverted: E = E ± (S ± ); E + is an 
increasing function, while E~ is a decreasing one. 

Each A* implies aperiodic family of trajectories A^ := A* + (2%j,Q) and corresponding quasimodes (£ ,± , Yj) 

with \ffj(x,h) = exp(27Zijx/h)\i/ ± (x,h) = T x J \j/ ± (x,h). To construct a Bloch quasimode we use an ansatz similar 
to that we use in the previous subsection, 

^(xXh) = £ cfyf{x,h) = £ cf r->±M). 

jez jeZ 

Therefore, Cj +1 — C^e 271 '^ 1 (C ± is a normalizing constant), and S ± (E ± ) — h(n^fk.2), n G Z. Denote 
E^(k2,h) : — E ± (S ± — h(n T^)); these functions can be viewed as semiclassical dispersion relations. Clearly, 
the functions are decreasing functions of k%, while E± are increasing ones, therefore, in some critical points k\ 
one has E* := E~(k%,h) = E^(kX,h), as illustrated in figure^. The corresponding points E* are usually treated 
as approximations of gaps, more precisely, one expects that in o(h°°) -neighborhood of each such value there is a 
gap, whose width is also o(h°°) [5]. The asymptotics of the true dispersion relations can be combined from pieces 
of and E m (figure|2b). 

As we see, we have o(h°°) -narrow bands in the first case and o(h°°) -narrow gaps in the second case. Our aim is to 
describe the transition between two these two extremal cases, which correspond to the semiclassical asymptotics 
near separatrices. To obtain at least the first non-trivial term in the asymptotics of the width of <9(/i°°)-smail 
bands and gaps one should take into account the interaction between neighboring cells, which is an extremely 
different problem [17]. Such a calculation involves topological characteristics of the Hamiltonian, which results in 
a description of the quantized Hall conductance [13]. It is interesting to emphasize that the resulting topological 
numbers come from index-like characteristics of a certain path on the Reeb graph [5]; tjis provides an additional 
path index-like interpretation of the Chern classes, related to their their description through the cycle related to 
edges states on the Riemann surface of the Bloch functions [19]. 

3. Singular Bohr-Sommerfeld rules and their modification for the periodic problem 

In this section, we give a short description of the semiclassical asymptotics near separatrices; this technique was 
developed in [14]. 

Let H(p,x) be an arbitrary classical Hamiltonian (not necessary a Harper-like one) with non-degenerate critical 
points and H/, be the corresponding quantum Hamiltonian. Let £ be a critical value of H in the sense that the 
level set A := {(p,x) G R 2 : H(p,x) = E} contains a saddle point of H; our aim is to study the asymptotics of the 
spectrum of Hh in the interval [E —Ah,E+Ah],A > 0. We assume that A is a compact connected set. 

3.1. Preliminary constructions. Near the critical energy it is more convenient to use the scaled energy, i. e. we 
are going to solve the equation 

(H h -E -Xh)^ = o(h m ), (7) 
where A is a new spectral value to be found; we consider the situation when X runs through some finite interval 
[—A, A], A > 0. The conditions which guarantee the existence of such solutions are called the singular Bohr- 
Sommerfeld rules. Due to the non-degeneracy of the saddle points, A is a tetravalent graph embedded into the 
plane R 2 , and the saddle points are point of branching. The edges of the graph are smooth curves; each of these 
curves, v, delivers a part of solution by means of the usual WKB-asymptotics [18]; we denote this function by 
y/ Y . Our aim is to glue these contributions together near saddle points (vertices of the separatrix graph) in order 
to obtain a requested solution. Let us introduce a enumeration of the edges of the separatrix near each vertex 
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in the following manner: the direct cyclic order is Vi,V4,V2,V3, and the quadrant formed by the edges v>i and V3 
is pointed to the top, see figure [5] Near each vertex V we try to construct the requested solution in the form 
\P = xj \ff vi Yv2 ®xJYv 3 ®xJy V3 , where xj e C, j = 1,2,3,4. (Sometimes we omit the superscript V and 
write simply Xj.) Therefore, an edge can have different numbers near different vertices. To any vertex V of the 
graph we assign a so-called semiclassical invariant £ V which is a formal power series in h, 

e Y (l,h) = YeJ(l)hJ, where e v (A) = ±^ A H" = (""" H J X 

Pa 1 ' ^|det//"(V)|' \Hxp H xx 

and the sign ± coincides with the sign of H in the quadrants formed by the edges 1 and 3. 

Each cycle 7 on the graph A will be accompanied by the following three characteristics (see [14] for details): 

• Principal action Ay, 

Ay=<j>p dx, 

• Renormalized time Iy. For cycles 7 crossing critical points with corners, we put 

Iy = v.p.<fdt:=Y[ J dt, c,-ey, ci = c n+ i, 

where the points Cj are chosen in such a way that each of the pieces (c/,C; + i) C 7 contains exactly one 
corner V/, and the integrals are calculated as 

dt:= lim / dt+ / dt + sAog 

a,b->V, \J C J h 



dpAdx 

fiab 



where Sj :— =FA/^/|det//"(V/)|, and R a \, is a parallelogram spanned by the points a, Vy, and b, and the 
sign =F is "-" if the direction of integration corresponds to dt and "+" otherwise. By dt we denote the 
Hamiltonian time, dt (sgrad//) = 1. This definition of Iy is then extended by additivity to all cycles (not 
necessary with corners). 

• Maslov index niy. If none of the tangent to 7 at the corners is parallel to the p-axis, the index my is 
calculated as the sum of the indices of all the turning points in which 7 is smooth. If there are tangents 
parallel to the p-axis, one can destroy them by small variation; this must be reflected in the enumeration 
of edges near all the vertices. 

3.2. Singular Bohr-Sommerfeld rules. The singular Bohr-Sommerfeld rules, which finally will result in condi- 
tions for A, come from the following procedure. 

Near each saddle point V there exist a canonical transformation %{q.,y) = (p,x), %{0) = V, such that 
H(p,x) = W(qy)qy with a certain smooth function W. This implies an elliptic Fourier integral operator U and 
a pseudodifferential operator W, elliptic at the origin, such that HU = UW(yq — he^). Using this representation 
one can construct a microlocal basis \ffj, j = 1,2,3,4, of semiclassical solutions near V. More precisely, we put 



«p lW =jB v ^l e '' el °sM (p 2 (y) -B v ^^ e '- £l °sl>'l 




«p 3 ( y ) = B v l= / J*/" e -^\'\dt : (p 4 (y) = B V ^= / e^" e'^dt, 




(8) 



where Y is the Heaviside function, and then Yj = U (pj is a basis element corresponding to the edge j, j= 1 , 2, 3,4. 
Here B v is a normalizing constant. 
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A linear combinationxi 1/1 (Bx 2 V2 ©*3V3 0*4 y/4, xj 6 C, j = 1,2,3,4, which extends the WKB-solution to the 
saddle point, defines a function near V iff x/s satisfy the linear system 

% )=<?(. L ie 7 n )( Xl ), *:= , 1 ^^H****, £ = £ v (9) 
K x A ) \ie £% 1 J \x 2 J ^/ l+e -2K£ 

Clearly, jt/'s corresponding to different vertices must be connected with each other in a certain sense. To describe 
this correspondence we cut several edges in order to get a maximal tree on the separatrix. Now consider an 
arbitrary edge v between two vertices V' and V". Denote the corresponding x-coefficients by x J j and x" respectively, 
j = 1,2,3,4, and assume that v has index j with respect to V' and index k with respect to V". We put 

x'j = x'[, if v is not cut, x 1 : = hoi yx", if v is cut, (10) 

where holy is the so-called holonomy of the cycle y formed by the edge e and the edges of the maximal tree (this 
cycle is unique). For the holonomy holds the estimate |holy| = 1, Argholy = Ay/h + Xly + nm-y/2 + 0(h). Note 
that the first equality in JlOi fix the constants B v uniquely (up to a common multiplicator). 

The equations (|9jl and QOj l written for all the vertices and all the cycles respectively compose a linear sys- 
tem for the coefficients xj, j — 1,2,3,4. The condition for the existence of non-trivial solutions (non-vanishing 
determinant) is called the singular Bohr-Sommerfeld rules. 

3.3. Periodic problem. Our aim is to apply the singular Bohr-Sommerfeld rules to the Harper-like operators. 
The problem is that the level set of E is always periodic and, respectively, unbounded. To apply this technique to 
the problem in question we take into account the Bloch conditions already at the stage of the construction of the 
solution (in contrast to the smooth case), in other words, we are going to apply the procedure described above for 
constructing the quasimode ( v t , (x,k,/i),£'(k,/i)) such that ^ satisfies ([Q and 0. Clearly, this condition can be 
rewritten as the set of equalities B v + 2jtl = B y , and xj +2nl = e^xj, j = 1,2,3,4, for all vertices V and 1 e Z 2 , 
where B v are constants from This clarifies the meaning of the maximal tree in the periodic case: Instead of 
constructing a maximal tree in the whole plane it is sufficient to find a maximal tree on the reduced phase space 
T 2 X ; all the cycles will be viewed then as cycles on the torus. 

4. Transition between finite and infinite motions 

In this section, we consider spectral regions corresponding to the transition between closed and open trajectories. 
In this case the separatrices are non-compact in one direction only, and we assume that they are directed along the 
x-axis. Clearly, one can proceed first exactly in the same way as in the case of open trajectories, which will result 
in the following ansatz for the requested solution: 

¥(*,k,/*) = c £ e iJkl T^ J y(x,k 2 ,h), 

where y/must be a quasimode associated with the separatrix, (H^—E — Xh)\j/= o(h°°), and satisfying 72 y = e lkl \y. 
It is clear that the corresponding semiclassical dispersion relations will not depend on £1 . 

4. 1 . Two edges of infinite motion and one edge of finite motion. In this subsection, we consider probably the 
simplest structure of the separatrix. More precisely, we assume that the energy level H = E on the reduced phase 
space T 2 px contains exactly one critical point, to be denoted by V. The corresponding separatrix on the plane M 2 T 
has then the shape showed in figure|4] In terms of the Reeb graph this situation means that we consider a transition 
between two edges of infinite motions and one edge of finite motion. Near the corresponding branching point the 
Reeb graph has a T-like shape. (We assume that the closed trajectories lie under the critical energy level.) Denote 
the semiclassical invariant of V by e(A,/z). Clearly, 



e(l,h)=l/w + 0(h), w:= >/|det#"(V)|. (11) 

The separatrix on T 2 px is a graph with one vertex and two edges, j\ and y>i, which are cycles on the torus, see 

figure^ In order to obtain a maximal tree one has to cut both edges. Denote hoi y, =: e ia > , j — 1,2, where 

A: . nnij . . f f 

aj = -± +Alj+ — - +0(h), Aj:=f pdx, Ij = v.p. at, m 7 =indy-, /*=1,2. 
h 2 J Yj ■ J 7j 



(a) 



(b) 



Figure 4. Transition between one edge of finite motion and two edges of infinite motion: (a) 
Structure of the separatrix, (b) The corresponding part of the Reeb graph 



To simplify the notation we put Xj :— xj, yj :— xj" rv "' z '"'', j = 1,2,3,4. The singular Bohr-Sommerfeld rules and 



V+(0,2tt) 

the Bloch conditions lead us to the following equations: 

Bloch-periodicity conditions: yj 

' x i = gx l +ie~ £K gx 2 
X4 — ie~ £n Sx\ +$X2, 



e^xj, j = 1,2,3,4, 



Matching at V: 



Holonomy equations: 



'y 3 = e icll xi, for 7!, 
x 4 = e' a2 y 2l for 72. 



where 



S : 



1 



O ;'argr(^+;'e)+!£log/! 



The condition of existence of non-zero solutions is equivalent then to the equation 



w 1 . \ , , «1 +«2 

cos( argr(- + is) + e\ogn — 



1 



-2ne 



cos ( k 2 — 



0,2 - tt\ 



The spectral parameter A enters this equation through e, ot\ and 0,2- One has obviously 



_/l \ CC\ ~t~ OC2 1 

aigL[- + ie) + e\ogh — = ±arccos^= — 

2 2 y 1 



,—2%z 



cos ( IC2 



Cbi-CL\ 



■2%n, n € Z. 



To find approximate solutions we represent A as A = Ao + jit with fl = o(l); in this way we can solve the equation 
near each Ao. Taking into account Eq. il It one can write 

(\ iXo\ Ao U (Ai+A 2 ) /1+/2, n(m l +m 2 ) 
argr(- + — +— log + - log A — Ao 

V2 w J w w In 2 2 

± arccos — cos (k2 ^ — - — ^ — -Aq") + 2nn + 0(A) + <9(/ilog/i), n e Z 



yl + e - 2K *o/» 
(one can show easily that m\ — 1112 = 0), or, finally 

w 



2h 



H = H„(k 2 ,Xo,h) 



log/1 



cos (k2 — A(Xo,h)) 
N(Xo,h) ± arccos v )= + Inn 



log/i 



N(Xo,h) ± arccos 



■J 1 -|-g-2fl;Ao/w> 
cos fe - A(Ao,A) 

y 7 J _|_ e -2ffAo/w 



1 + 



27T« 



log/i 



■0(A) 



(12) 



where 



1 \ A i+ A 2 . A (/i +I 2 ) + 7t(mi +m 2 ) Aolog/; 
N(Xo,h) = 1 argT 



2/7 



A(Ao,/z) = 



2 



1 | 'Ap 

2 w 



Ao 



/2-/1 



2/z 2 

and « are integers such that the expression in the curly brackets in dl2l is o(\ogh). Returning to the original spectral 
parameter we obtain a series of semiclassical dispersion relations: in a o(h) -neighborhood of E + XqIi they take the 




(a) (b) 

Figure 5 . Transition between two edges of finite motion and two edges of infinite motion: (a) 
Structure of the separatrix, (b) The corresponding part of the Reeb graph. 



form E„(ki,k2,ha,h) = E + X^h + hjj.^(k2,Xo,h). The expression obtained can be used for estimating the band- 
and the gapwidth. More precisely, the bands and the gaps in a o(h) -neighborhood of E + XqIi have the width 

VI Vlogrt/ 



2wh 

B(Ao,h) = - — -arcsin 
| log/* I 

G(Xo,h) = — — arccos 



1 



= + °(A) 
w V log h I 



(14) 



| log /z | y/l+e- 2 "*o/ 

respectively. In particular, for Xq = the bands and the gaps have approximately the same width nwh/(2\logh\). 

4.2. Two edges of infinite motion and two edges of finite motion. We consider now another structure of the 
separatrix. Let us assume that the energy level H = E on the reduced phase space is a connected set containing 
two critical points, which we denote by V and V. The separatrix has the shape sketched in figure [5] The Reeb 
graph has near the corresponding branching point a X-like shape, where the two upper edges correspond to open 
trajectories. 

Denote the semiclassical invariants of V and V by e and e respectively, 



e(l,h)=l/w + 0(h), 



■ k/w+0{h), w: 



det//"(V)| 



w := 



detff"(V)| 



and put 



£': = 



exp 



/argT( \ +ie ) +/elog/i, 



exp 



i'argT^ + + /elog/z, 



-2nl 



(15) 



(16) 



In order to obtain a maximal tree on the reduced phase space we cut all the edges but j\, then one get three cycles: 
fi = Yi + 72, 72 = 74 - 7i> and 73 = 7i + 7?; put CCj : = argholy,-, j = 1,2,3. 



Denote Xj := xj, jj := Xj, z 



_Y+{0,2n)^ j _ 1 2^3,4, then the quantization conditions are: 



Matching conditions at V: 



; ' ~J ■ "j 

Bloch-periodicity conditions: ; 

X3 = §X\ + ie 
£n Sx 



— e ikl x ■ 



A'4 = ie 



+ 3X2, 



at V: 



j= 1,2,3,4, 
'y 3 = iyi + ie-' ex iy 2 , 



y 4 = ie 



Holonomy equations: y>2 = X4, yj = e l X\ 



Z3 = e *yi, y 4 = e 3 z 2 - 



The condition of the existence of non-trivial solutions leads to the equation 



cos j argT^i + ie 



argr^ + /e 



(e + e)log/i 

cos ( k2 



CC1 + CC2 + 0C3 



'{l+e- £7Z ){l+e- £71 
The parameter A enters this equation through e, £, and aj, j = 1,2,3. 



-fe+S)n a l — a 2 — «3 
-e { ^ ' COS 



(17) 
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To obtain a more clear picture, let us introduce holonomies of the edges as solutions of the following equalities: 
OEi = ft + ft, OC2 = ft — ft, aj, — ft + ft, arghol(72 + 74) = ft + ft; these solutions can be represented as 

6, „ Km; , , f 
P j = -±+XJ j + ^ L + 0{h), Bj = J^pdx, mj =mdyy, 7 = 1,2,3,4, 

then illl takes the form 

ft+ft + ft + ft 



r g r(i + /e) +argr(i + /e) + (e + £)logft 



2 



COS (* 2 - fe&HMl - e^+^COS (^+fe)~(ft+ft) 

= ±arccos ^ I h27T«, neZ. (18) 

The quantities ft can be viewed as "weights" of the separatrice edges, then Eq. Jl 81 shows how the relationship 
between them influences the dispersion relations. 

Like in the previous case we represent A as Ao + jU, 11 = o(l), and solve (approximately) the equation for 11. 
The solutions take the form 



(19) 



where 



1 WW 

ll = Li±(k 2 ,lo,h) = - — — — — - x <^ N(lo,h) 
[w + w) log ft y 

cos(^~A 2 (Ao,/i))- e -( 1 /»'+ 1 /^cosA 1 (A ,/i) 1 / , 1 x\ 
iarccos , , h27Tn > x 1 + - — - +0(h), 

^( 1+e -2^ /w)( 1+e -2^/*) J V Vlog/lV V h 

»r/o ,\ 61+62 + 63+64 , , 7i+72 + /3 + -/4 , w(mi +m 2 + m3+m 4 ) 
:= ^ + A + 

- ar g r(i + ^) - ar g r(i + ^) - Ao(- + -) log A, 

(61 +6 2 )- (63+64) (J l+ J 2 )-(J 3+ J 4 ) 
Ai(/o,AJ .= — hAo ™ ' 

A ... (62+64) -(6l +63) , . (Jl+J4)-(Jl+h) 

A 2 (Xo,n) := ^ hAo 

and n 6 Z, « = [iV(Ao,ft)] + o(logft). The dispersion relations are E{k\,k2,h) = £ + /1A0 +hn^(k2,Xo,h). In 
contrast to the previous situation, band and gaps have, generally speaking, different width. More precisely, near 
the point E + XqIi one has bands having the width 



. „ , wwh 
6(Ao,/i): 



(w + w)|logft| 



arcsin ■ 



1 +e -(l/"'+l/w)Ao7r cos Ai(Aq,/») 
v/(l + e - 2 ' I V»')(l + e -27tV*) 

l_e-(i/»'+i/w)Aoff cos Ai(Ao,A) 
+ arcsin — — 
a/(1 + e -23tVM-)(l + e -2^/*) 

and two groups of gaps having the width 



/ ft 

CM 

V log ft 



(20) 



n 2wwh l+e-a/v^+i/^cos Al (Ao,/0 , / ft \ 

Gi(/io,/i) = - r-: rarccos — |-o , (21) 

(w + vv)|logft| y/(i _)- e -2wAo/x')(l + e -2wAo/*) Vlogft/ 



2^ i_ e -(i/>^+i/»)^ cos Al (Ao ,ft) / ft \ 

G2(M).h) = rr -arccos — - +0 , (22) 

(w + w)|logft| +e -2^Ao/»')(l + e -2^V*) Vlogft/ 



and they come in groups . . . ,6, G\ ,6, G2, . . . , The formulas clearly show the fast decaying of the bandwidth for 
negative Ao and of the gapwidth for positive Ao- From the other side, the ratio bandwidth/gapwidth depends 
crucially on the relationship between Bj, Jj, Ao, and ft. 

5. Transition between topologically different finite motions (degenerate case) 

In this section we consider a situation when all the smooth trajectories of the classical Hamiltonian are closed. 
This situation takes place if, for example, H is invariant under linear transformation having no real eigenvectors 
(rotation by n/2, for example). More precisely, we assume that the level setH = E on the reduced space contains 
two critical points, which we denote by V and V, and the corresponding separatrix in the plane is non-compact 
in all directions. This situation is illustrated in figure[6] We can expect that the semiclassical dispersion relations 
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(a) (b) 

Figure 6. Transition between two edges of finite motion: (a) Structure of the separatrix, (b) 
The corresponding part of the Reeb graph. 



in this case depend on both quasimomenta k\ and k%. Denote by e and e the semiclassical invariants of V and V 
respectively, and use the notation of dl6> . Clearly, 



-k/w + 0{h), w 



det//"(V)| 



dettf"(V) 



An essential point is that the main terms of e and £ have opposite signs. 

We fix a maximal tree by cutting all the edges but J\, then three cycles appear: fl = 74 + 7i, 72 = 73 — 7l> an d 



V+(0,2?r) 



72 = 72 + 7i with holonomies e ,a >, j = 1,2,3, respectively. Put xj := jcJ , y j : = Xj , Zj '■= x- 
j = 1,2,3, 4, then we come to the following set of equalities: 

Bloch-periodicity conditions: w j — e' kl yj, Zj = e Xj, j = 1,2,3,4, 



, Wj 



V+(27T,0) 



Matching conditions at V: 



' x 3 = Sx\ + ie~ en £x2, ~ { y3 = &y\ + ie~ £7t iy 2 , 



Holonomy equations: X4 = y$, W2 — e l x\ 



xa, = ie £n §x\ + £x2, 



at V: 



y 4 = ie- £K £ yi +£y 2 , 



z 3 = e' a2 w 4 , 



After some algebra one arrives at a 2 x 2 linear system, 

g — £> e i{ai+ai-h) ie~ £7l g — i e -B% g e i{ki+a l +a i ) 



e ia 'z 2 . 



= 0. 



The condition for the last system to have non-trivial solutions comes from the vanishing of its determinant and has 
the form 



cos 



argrQ + ie) -argr(--Me) +(e-g)logA- 



«1 + Oil + «3 



-2£K 



)(1 



-2ex\ 



«-( e+e )*cos(* 



- cos ( k 2 



0,1+0,2 — 03 



(23) 



Introducing again the weights j3, of the edges y,-, j — 1 , 2, 3,4, by the rule 

ai = j8 4 + j8i, a 2 = j3 3 -j8i, a 3 = j3 2 -j3i, arghol(y 2 + y 3 ) = /3 2 + /3 3 , 
B j . Ttm j 



Pj = -^+kJj + ^ + 0(h), Bj 



IS 



pdx, mj = indYj, 7 = 1,2,3,4. 
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we rewrite d23l as 



COS 



argr(- + /e) -argr(- + /e) + (e-£)log/z- 



A 



(MP±LdMM\ J. cos (k 2 - (ft+foHPi+ft) 



2 -r^i^z 2 



where 



(24) 

'l +e- 2£7I )(l +e- 2lK ) 
Representing again X = Ao + jtt, /i = o(l), we come to the following expression for /x: 

i WW 

M = M„ (^i,^2,Ao,/j) = t — — - — - 
(w + wj logn 

f e -^cos(^-A 1 (Ao,/ i ))+ e -^cos(^ 2 -A 2 (A ,/ ? )) 1 
x |^,*) ±arccos ^^7^^ ^ 

iV(Ao,ft) = + Ao + 

- (1 + i)A log/, - argrfl + ^) + argrfl - ^ 

ww \2 w J \2 w 

(B l+ B 4 )-(B 2 +B 3 ) (/ 1+ y 4 )-(7 2+ 7 3 ) 

Mh,h) = Yh + *° 2 ' 

( j B 3 +fi4)-( J Bl+fi 2 ) , . (-/ 3 +/4)-(/l+/2) 

A 2 (A ,n) = — hAo . 

2n 2 

The bandswidth B(Ao,/z) and the gapwidth G(Xo,h) in a o(h) -neighborhood of the point E + XqIi admit a simple 
estimate 

B(Xo,h) — r arcsin— ■ — = +0 , (26) 

(w + w)|l0gA| ^/(i +e -27tAoA")(l +e -27rAo/»f') Vlog/l/ 

2ww/! g-»V*4- e -' t V* / h \ 

G(AQ.h) = r-; rarccos — = +0 , (27) 

(w + w)|logrt| _|_ e -2^Ao/ ,v )(l +e -2 t^o/*) Vlog/i/ 

so that the bands clearly show an exponential decay with respect to |Ao|. 

6. Discussion 

In this section, we discuss in greater detail the influence of characteristics of the Hamiltonian H on the dispersion 
relation. 

In the case of subsection !4. li the picture is quite simple. The semiclassical dispersion relations depend on one 
of the quasimomenta only, and the extremum points (with respect to this quasimomentum fe 2 ) is determined by 
the difference between the upper and the lower part of the separatrix, j\ and y 2 ; the width of the bands and the 
gaps, which is given by Eqs. Jl 3i and (I14> respectively, depends on the determinant of the second derivatives at 
the critical point; these formulas present a generalization of a similar estimate for the periodic Sturm-Liouville 
problem (see, for example, [14, Sec. 10]). 

The example considered in subsection 14.21 and figure [5] differs from the previous one. The position of the 
extremum point with respect to fc 2 is, like in the previous case, determined by the relationship between the upper 
72 + 74 an d the lower j\ + 73 parts of the separatrix. But the band- and the gapwidths, as can be seen from J20I . 
J21L and I22L crucially depend on Ai; the quantity A\ can be interpreted as a "difference" between the cycles 
7i + 72 an d 73 + 74- F° r example, if the areas of these two cycles do not coincide, the ratio bandwidth/gapwidth has 
no limit for h — > 0. Examples with more complicated separatrices will show a more curious picture. 

The example of section|5]shows a regular behavior of the band and the gaps; the main term in the asymptotic of 
their width, Eqs. (I26> and (I27> . like in the first example, is determined by the second derivatives at the critical points 
only. From the other side, this example is suitable for discussing the form of the energy bands. The dispersion 
relations lying in a neighborhood of E + XqIi have maxima at the points 

k = W^*):=(2»x{^^ + I}- W> 2 W x{^^ + !}-«) (28) 
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Ao>0 Ao = Ao<0 

Figure 7 . Levels curves of the dispersion relations lying near E + Aq/z 



(here {x} denotes the fractional part of x). In the generic situation, when all the holonomies j3 ; - are different, these 
points depend crucially on Ao, Bj, Ij, and h, so that a small variation of them can change k max significantly. Near 
the critical energy E (i.e. for Ao = 0), the quasimomenta have "equal rights", i. e. the coefficients before cos-terms 
in (I25> are approximately equal to 1 . In the non-symmetric case, when w ^ w, the situation changes for non-zero 
Ao- So, if w > w, the quasimomentum k\ dominates if Ao > and ki dominates for Ao < 0, and, at the same time 
the maxima of the dispersion relations moves according to ( 128k see illustration in figure Q This gives a rather 
rough (but at the same time generic) impression about the dispersion relation structure near the critical point. 

Appendix. The periodic Landau Hamiltonian with a strong magnetic field and Harper-like 

operators 

In this appendix, we briefly describe the relationship between the periodic Landau operator and Harper-like 
operators as it was established and studied in [4, 1 1, 20]. 
The periodic Landau Hamiltonian has the form 



L:= 



-ih 



dx\ 



-X2 



■ev(xi,x 2 ), 



f/o(v / -2/A 3 ,)v(yi,y 2 ), 
One can show that there exists a 



2 dx2 

where visa two-periodic function with periods a and b. If £ = 0, the spectrum of L consists of infinitely degenerate 
eigenvalues /„ = (n+ 1 /2)h, n G Z+, called Landau levels. The presence of non-zero £ leads to a broadening of 
these numbers into a certain sets called Landau bands. We are going to show, under assumption that both h and e 
are small, that the broadening of each Landau level is described by a certain Harper-like operator. 

The corresponding to L classical Hamiltonian is L(p,x,e) = (p\ +X2) 2 /2 + p\j2 + Ev{xi,X2)- If £ = 0, then 
L defines an integrable system whose trajectories on the (x\ .X2)-plane are cyclotron orbits. For non-zero £ the 
Hamiltonian is non-integrable, but, for small £, one can interpret the classical dynamics as a cyclotron motion 
around a guiding center. We introduce new canonical coordinates connected the motion of the center: p \ = —y2, 
P2 = —q, x\ = q + yi, X2 = p+y2, considering (p,yi) as generalized momenta and (q,y 2 ) as generalized positions 
(p, q describe the motion around the center with the coordinates y\ ,y 2 ), then H takes the form L — (p 2 + q 2 )/2 + 
Bv(q+yi,p + y 2 ). Introduce the averaged Hamiltonian 

L(I,y u y 2 ,e) =1+ — Iv(V2Ism<p+yi,V2lcos<p+y2)d(p = , 

where Jq is the Bessel function of order zero and A, = d 2 /dy\ + d 2 jdy\ 
canonical change of variables (p,yi,q,y2) = (P,yi,q,y2) + O(e), periodic in y\ ,y 2 with periods a and b, such that 
L = L{ y (p 2 +q 2 )/2,yi 7 y2,e) +0(s 2 ). The averaging procedure can be iterated, so that one constructs an averaged 
Hamiltonian Jzf = J£(J,Y\,Y2,e) and a canonical change of variables (p,yi,q,y2) = (P,Yi,Q,Y2) + 0(e), both 
periodic in yi,y2, such that L = JC((P 2 + Q 2 )/2,Yi,Y2,e) + 0(e°°). Therefore, neglecting the last term and using 
the canonicity of all the transformations one reduces the spectral problem for L to that for J?f obtained from Jzf by 
the Weyl quantization (P = -ihd/dQ, Y\ = -ihd/dY 2 ): 

&¥{Q,Y 2 )=EV(Q,Y 2 ). (29) 

Clearly, commutes with the harmonic oscillator / := (—h 2 d 2 /dQ 2 + Q 2 ) /2, which means that the eigenfunction 
*¥ in d29l can be represented as ^(Q^) — \j/ n (Q)<P(Y2), where y/„ is an eigenfunction of / with the eigenvalue /„, 
and <t> must be an eigenfunction of the operator L„ obtained by quantizing the classical Hamiltonian L n (Y\,Y 2 ) — 
J*?(I n ,Y{,Y2,£) considering Y\ as a momentum and Y 2 as a position. All these Hamiltonians L n are periodic in Y\ 
and Y2, therefore, after a linear symplectic transformation L„ becomes a certain Harper-like operator which can be 
treated as a Hamiltonian describing the broadening of Landau level under the presence of the electric potential v. 

An essential point in our considerations is the dependence of L n on n: one has L n (Y\,Y2) = I n + 
£Jo{y/— 2/„Ay)v(7i ,Y 2 ) + 0(e 2 ). In general, the topology of the trajectories of L n depends on n, and a given 
potential v can produce a number of operators L n defining the structures of different Landau bands. Considering a 
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Figure 8. Level curves in the elementary cell [0,2k) x [0, n): (a) of the potential v, (b) of the 
averaged Hamiltonian for / ss 0.5, (c) of the averaged Hamiltonian for 1.5 

simple example v = cos 2 (xi/2)cos 2 X2 one arrives at Jzf(I,Yi,Y2,£) = I+s(l +Jq(V2I) cosxi +7o(v / 87)cos2x2 + 
-A)(v'10/)cosxi cos 2x 2 ))/4 + 0(e 2 ). The level curves of the potential v and of the averaged Hamiltonian are 
sketched in figure|8l obviously, they are different, which implies differences in the structure of the corresponding 
Landau bands. 
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